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ABSTRACT 


This  thesis  examines  the  problem  of  predicting  distortion 
induced  stall  in  an  axial  flow  compressor.   Three  instantane- 
ous distortion  indices  are  used  on  the  data  presented.   The 
data  consist  of  eleven  stall  events  run  at  2.5  Mach  by  NASA 
Lewis  on  a  J85-GE-13  turbojet  engine.   The  three  indices 
used  include  two  indices  based  on  vorticity  and  developed  at 
the  Naval  Postgraduate  School.   The  third  index  was  developed 
by  Pratt  and  Whitney  and  is  based  on  total  pressure  distri- 
butions.  The  three  indices  are  compared  by  visual  means  us- 
ing computer  plots  and  are  also  compared  by  using  a  simple 
statistical  approach.   The  computer  program  used  was  first 
developed  by  Shoemaker  in  his  Master's  thesis.   The  program 
is  expanded  in  this  thesis  to  compute  and  compare  the  three 
indices. 
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I.   INTRODUCTION 

This  paper  addresses  the  problem  of  distortion-induced 
or  drift  stall  in  an  axial  flow  compressor.   The  problem  is 
a  current  one  in  industry  and  has  not  been  wholly  resolved. 
The  present  method  of  dealing  with  the  problem  of  distortion 
is  to  operate  axial  flow  compressors  with  a  sufficient  stall 
margin  to  avoid  the  problem  altogether,  but  this  is  an  expen- 
sive procedure.   Best  engine  efficiencies  are  obtained  by 
operating  close  to  the  stall  or  surge  line.   By  providing  a 
large  stall  margin,  we  are  in  effect  sacrificing  performance 
via  reduced  efficiencies.   The  problem  is  especially  germane 
to  military  aircraft.   The  normal  operating  conditions  of 
military  aircraft  include  hard  maneuvering,  high  "g"  flight, 
steam  ingestion  from  aircraft  carrier  launches,  and  exhaust 
gas  ingestion  from  forward  firing  weapons.   These  operating 
regimes  cause  a  pressure,  temperature,  and  gaseous  molecular 
weight  distortion  of  the  flow  entering  the  compressor  section 
of  the  engine. 

In  order  to  deal  with  distortion,  it  is  first  necessary 
to  find  or  define  an  adequate  measure  of  it.   A  distortion 
index  based  on  steady-state  or  long-time  averaged  values  of 
distortion  has  been  attempted.   This  method  has  met  with  some 
degree  of  success,  but  it  is  apparent  that  peak  instantane- 
ous distortions  of  sufficient  magnitude  can  induce  stall  when 
the  time  averaged  distortion  factor  remains  relatively  low. 
The  current  approach  and  the  approach  of  this  paper  is  to 
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examine  distortion  as  an  unsteady  phenomenon  and  to  measure 
distortion  over  time  intervals  comparable  with  the  character- 
istic times  of  an  axial  compressor. 

A  computer  program  was  first  developed  by  Ens.  Shoemaker 
[ref  3]  utilizing  an  existing  interpolation  scheme  developed 
at  the  Lewis  Research  Center  of  NASA.   The  program  calculates 
the  total  pressure  distribution  just  prior  to  the  compressor 
section  of  a  J85-GE-13  turbojet  engine  as  a  function  of  time. 
The  program  also  utilizes  the  theory  of  Fuhs ,  Farmer,  and 
Iverson  [ref  1,  ref  2]  to  compute  the  circumferential  and 
radial  components  of  vorticity  across  the  face  of  the  com- 
pressor section  as  a  function  of  time.   The  computer  program 
was  extended  for  this  paper  to  calculate  and  statistically 
compare  two  existing  distortion  indices  plus  an  additional 
index  developed  by  Downs  and  studied  in  this  thesis  [ref  4]. 
This  procedure  was  applied  to  eleven  stall  events  recorded 
and  analyzed  by  the  NASA  Lewis  Research  Center  [ref  5]. 
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II.   NATURE  OF  THE  PROBLEM 

A.   COMPRESSOR  SURGE 

The  line  marked  "surge"  in  figure  1  denotes  the  locus  of 
unstable  operation  points  for  an  axial  compressor.   The 
cause  of  compressor  surge  may  be  visualized  in  the  following 
manner.   At  any  given  reference  RPM  (N//6")  and  blade  radius, 
a  rotor  speed  U  is  specified.   Consider  the  velocity  diagram 
in  figure  2.   A  reduction  in  referred  mass  flow  rate  will 
cause  a  reduction  in  the  axial  component  of  velocity.   Since 
the  rotor  speed  is  fixed,  an  increase  in  the  relative  flow 
angle  3  will  occur.   Present  compressor  design  incorporates 
a  fixed  rotor  geometry  and  thus  the  relative  flow  angle  may 
only  increase  to  a  maximum  value  which  is  dependent  on  blade 
design . 

An  insight  into  the  unstable  nature  of  the  surge  line 
may  be  made  by  comparing  the  characteristics  of  the  compres- 
sor to  the  characteristics  of  a  throttle.   A  throttle  is  set 
to  produce  a  certain  back  or  receiver  pressure,  given  an  in- 
let pressure.   We  may  visualize  a  system  consisting  of  a 
compressor  producing  a  high  inlet  pressure  followed  by  a 
throttle  which  dumps  the  high  compressor  air  to  ambient  con- 
ditions (fig  3).   A  plot  of  the  characteristics  of  the  hypo- 
thetical throttle  and  compressor  are  included  in  figure  4. 
The  compressor  is  represented  at  a  fixed  reference  RPM.   The 
dashed  lines  represent  fixed  throttle  settings.   Selecting  a 
particular  throttle  setting  determines  two  possible 
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compressor  points.   Consider  as  a  throttle  setting  line  B, 
with  possible  operating  points  P  and  Q.   First  consider  point 
P.   If  the  mass  flow  rate  is  perturbed  to  a  higher  mass  flow 
rate,  the  throttle  setting  will  demand  a  higher  pressure 
ratio,  whereas  the  compressor  produces  a  lower  pressure  ra- 
tio.  In  this  case  the  combination  of  compressor  and  throttle 
will  return  to  point  P.   Conversely  a  lower  flow  rate  will 
cause  the  throttle  to  demand  a  lower  pressure  ratio  whereas 
the  compressor  will  produce  a  higher  ratio.   Again  the  sys- 
tem will  be  forced  to  return  to  point  P.   Consider  now  point 
Q.   A  slight  perturbation  mass  flow  will  produce  problems. 
The  system  is  free  in  this  case  to  follow  any  perturbations 
and  the  system  is  unstable.   Neutral  stability  is  represent- 
ed by  point  R  where  the  slope  of  the  throttle  curve  matches 
the  slope  of  the  compressor  curve. 

B.   CLEAN  STALL  VERSUS  DISTORTION  INDUCED  STALL 

The  surge  line  in  figure  1  denotes  the  conditions  at 
which  stall  will  occur  if  the  inlet  velocity  profile  is  uni- 
form, producing  a  uniform  total  pressure  distribution.   In 
actuality,  the  velocity  profile  is  not  uniform  and  pressure 
perturbations  may  produce  a  local  condition  over  a  small  re- 
gion of  the  compressor  face  which  will  cause  stall  over  that 
region.   This  stalled  region  in  turn  produces  a  condition 
known  as  rotating  stall.   In  figure  5,  if  the  conditions  are 
such  that  blade  2  becomes  stalled,  the  effective  flow  area 
between  blades  2  and  3  is  reduced.   This  reduction  in  flow 
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area  causes  a  flow  diversion  which  tends  to  overload  and 
stall  blade  3.   Continuation  of  this  process  will  produce 
rotating  stall.   It  is  noteworthy  that  rotating  stall  propa- 
gates in  the  opposite  direction  from  rotor  velocity  at  about 
half  rotor  speed.   Rotating  stall  may  develop  into  surge  or 
fully  developed  stall. 

It  is  apparent  that  the  instantaneous  velocity  distribu- 
tions across  the  compressor  face  may  cause  stall  when  the 
steady-state  conditions  are  below  the  surge  line.   This  ve- 
locity distribution  in  turn  produces  a  total  pressure  distri- 
bution which  is  the  form  of  the  data  examined  in  this  thesis 
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III.   DISTORTION  INDICES 

A.   VORTICITY  VIA  CROCCO'S  THEOREM 

Non-dimensional  vorticity  components  in  both  the  circum- 
ferential and  radial  directions  of  a  cylindrical  coordinate 
system  may  be  calculated  using  the  theory  of  Fuhs ,  Farmer  and 
Iverson  [ref  1,  ref  2).   Starting  with  Crocco's  Theorem  and 
using  an  order  of  magnitude  analysis,  it  was  possible  to  re- 
duce the  expression  for  these  parameters  to: 

-1    3P1    M  3T1      .       1     9P'     M   3T* 


6    yP'U'  8r'    2  9r'  '    r    yP'r'U'  89     2r'  di 
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D.    PT      ,    W0D      ,    2r 
PT      6     2a         D 

For  a  detailed  explanation,  consult  reference  1.   The  calcu- 
lations are  examined  more  briefly  in  reference  2. 

Analytical  studies  and  data  analysis  have  been  performed 
in  an  attempt  to  relate  the  effects  of  total  temperature  and 
total  pressure  distortions  on  compressor  performance  [ref  7], 
The  vorticity  approach  is  well  suited  to  this  problem,  as 
both  total  temperature  and  total  pressure  may  be  taken  into 
account  in  the  calculations  of  vorticity.   It  would  be  bene- 
ficial to  examine  data  which  consist  of  both  temperature  and 
pressure  patterns  across  the  compressor  face.   Additional 
work  in  this  area  would  be  helpful.   The  data  examined  in 
this  thesis,  however,  did  not  present  a  significant  tempera- 
ture gradient  and  the  expressions  for  vorticity  could  be 
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reduced  to 


i 


-1   9P   .    •       1     3P 


6    yP'U'  9r'  '    r    vP'r'U'  9  6 
z  z 

The  task  of  this  thesis  was  to  determine  whether  an  index 
based  on  vorticity  would  provide  a  meaningful  prediction  of 
compressor  stall.   Two  such  indices  based  on  vorticity  are 
examined  in  this  thesis.   The  first  index  was  proposed  by 
Shoemaker  [ref  3]  and  is  based  on  circumferential  vorticity. 
A  second  index  proposed  by  Downs  [ref  4]  and  based  on  both 
radial  vorticity  and  total  pressure  is  examined.   In  order 
to  provide  a  means  of  comparison,  a  third  index  based  on  to- 
tal pressure  and  taken  from  Zonars  [ref  8]  was  also  examined. 

B.   SHOEMAKER  INDEX 

James  E.  Shoemaker  approached  the  problem  of  finding  a 
vorticity  based  index  by  first  printing  out  maps  of  pressure, 
circumferential  vorticity  and  radial  vorticity.   In  order  to 
accomplish  this  goal  he  utilized  a  set  of  subroutines  devel- 
oped at  Wright  Patterson  Air  Force  Base,  and  altered  them  to 
print  out  the  maps  at  times  near  stall  and  at  times  far  from 
stall,  to  provide  a  comparison.   A  copy  of  his  computer  pro- 
gram is  included  at  the  end  of  his  Master's  thesis  [ref  3]. 
Using  this  procedure,  Shoemaker  found  a  vorticity  pattern 
that  frequently  occurred  within  twenty-five  milliseconds  of 
stall,  and  was  not  common  at  other  times.   The  pattern  was 
characterized  by  a  ring  of  large  negative  circumferential 
vorticity  at  the  midspan  of  the  compressor  face,  and  a  large 
positive  ring  at  the  tip.   He  developed  an  index  based  on 
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this  pattern  which  will  be  referred  to  as  the  Shoemaker  in- 
dex.  He  divided  the  compressor  face  into  three  concentric 
regions.   Region  one  at  r'  less  than  0.5  was  ignored.   Region 
II  was  established  as  that  area  between  radii  of  0.55  and 
0.70.   Region  III  was  that  area  between  radii  of  0.8  and 
0.95.   The  Shoemaker  index  then  summed  over  Region  II  but  set 
any  positive  circumferential  vorticity  equal  to  zero.   Like- 
wise Region  III  was  summed  where  any  negative  vorticity  was 
set  equal  to  zero.   The  instantaneous  Shoemaker  index  was 
set  equal  to  the  sum  of  the  absolute  value  of  the  two  summa- 
tions.  In  notation  form, 

Shoemaker  Index  (T)  =  r=  E  [w„]  -  r;  I  On] 

where    wfl  =  positive  values  of  circumferential  vorticity  in 
Region  III 
0)o  =  negative  values  of  circumferential  vorticity  in 
Region  II 
In  order  to  take  into  account  the  added  effect  of  the 
duration  of  the  pattern,  Shoemaker  set  his  final  index  equal 
to  the  product  of  three  consecutive  instantaneous  indices  de- 
scribed above.   This  latter  procedure  was  not  used  in  the 
present  thesis.   The  Shoemaker  index  as  used  herein  thus  re- 
fers to  the  instantaneous  index  and  is  an  indication  of  dis- 
tortion in  the  radial  direction. 

The  Shoemaker  index  has  one  real  advantage  when  the 
transport  time  from  the  measuring  probes  to  the  compressor 
stage  (where  stall  originates)  is  considered.   The  pattern 
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develops  at  a  time  before  stall  which  is  larger  than  the 
transport  time  from  the  measuring  probes  to  the  stall  region. 
This  could  be  very  useful  if  variable  stators  or  the  fuel 
control  were  available  to  avert  stall. 

C.   DOWNS  INDEX 

A  second  index  which  was  proposed  by  Downs  [ref  4]  was 
formulated  and  is  explored  in  this  thesis.   This  index  incor- 
porates both  total  pressure  and  radial  vorticity.   The  basic 
idea  was  to  consider  the  contribution  of  both  total  pressure 
and  vorticity  to  the  relative  flow  angle  3.   Increasing  3 
would  cause  the  rotor  blade  incidence  angle  to  increase  and 
would  tend  to  load  the  blade.   A  perturbation  plot  of  3  ver- 
sus 9,  where  the  deviation  of  3  from  the  reference  line  takes 
into  account  pressure  and  vorticity  effects,  is  shown  in  fig- 
ure 6.   The  corresponding  plot  of  total  pressure  and  to'  is 
also  shown.   The  reference  line  represents  uniform  flow.   It 
can  be  seen  that  although  go'  is  dependent  upon  the  pressure 
gradient,  the  combined  effect  of  pressure  and  to '  on  3  is  de- 
pendent on  both  the  gradient  and  relative  magnitude  of  pres- 
sure.  The  Downs  index  was  formulated  to  search  the  compressor 
face  for  the  largest  value  of  the  combined  effect  of  pressure 
plus  to',  over  any  60  degree  sector.   In  notation  form, 

P  -P 
Downs  Index  (t)  =  £  [~ — -  +  to'] 

P™        max  60 

where   P   =  average  total  pressure  over  entire  face  of 

compressor  at  time  t 
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The  Downs  index  is  primarily  an  indication  of  distortion  in 
the  circumferential  direction.   For  a  more  complete  analysis, 
consult  Downs  [ref  4]. 

D.   PRATT  AND  WHITNEY  K„„  DISTORTION  INDEX 

DM 

In  order  to  provide  a  means  of  comparison  between  vortic- 
ity  based  indices  and  an  index  calculated  directly  from  total 
pressures,  the  Pratt  and  Whitney  K  M  index  was  selected  from 
Zonars  [ref  8].   This  index  was  selected  primarily  because 
of  case  four,  which  was  one  of  the  eleven  cases  investigated. 
Distortion  indices  are  plotted  for  all  eleven  cases  in  Appen- 
dix A.   Neither  of  the  vorticity  indices  had  satisfactorily 
predicted  stall  in  case  four;  from  examining  the  pressure 

maps ,  K„,  appeared  to  be  an  appropriate  index.   The  K^.,  index 
DM  DM 

was  applied  to  case  four  and  to  all  other  cases.   This  index 
was  developed  by  Pratt  and  Whitney  for  the  TF-30  installed  in 
the  F-111A  aircraft  and  is  defined  as  follows: 

P     -  P 

i  v  r  Tmax Tmin,   _-  „ 

*  E  L     P~       Ji  ei  Ci 

v        _  Tave 

*DM  x  100 

i 

where    C  =  ratio  of  compressor  inlet  radius  to  ring  radius 

i  =  number  of  ring 

G   =  largest  continuous  arc  of  the  ring  over  which 

the  total  pressure  is  below  the  ring  average 

pressure 
PT    =  ring  maximum  total  pressure 

Pm    =  ring  average  total  pressure 
Tave  ^ 
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P_  .   =  ring  minimum  total  pressure 
Tmin       to  * 

It  can  be  seen  that  the  K^..  index  is  an  indication  of  dis- 

DM 

tortion  in  the  circumferential  direction  and  puts  more  empha- 
sis on  the  region  near  the  hub.   It  would  be  expected  that 
the  K  ..  index  and  the  Downs  index  would  be  comparable  since 
both  search  for  circumferential  distortion.   This  was  in 
fact  the  case  and  will  be  discussed  in  a  later  section. 
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IV.   EXPLANATION  OF  DATA 

The  NASA  Lewis  Research  Center  conducted  an  experimental 
wind-tunnel  investigation  to  determine  the  effects  of  time- 
varying  total-pressure  distortions  produced  in  a  supersonic 
inlet  on  a  J85-GE-13  turbojet  engine  [ref  5].   Selected  data 
taken  during  this  investigation  were  digitized  on  magnetic 
tape;  the  latter  is  the  source  of  data  used  in  this  thesis. 
The  following  explanation  of  data  is  largely  a  brief  summary 
of  the  explanation  presented  by  Burstadt  and  Calogeras  [ref 
5]  and  is  given  here  to  provide  an  overall  understanding  of 
the  data  analysis. 

The  inlet  was  an  axisymmetric  mixed-compression  type. 
The  compressor  was  an  eight-stage  axial  flow  compressor  with 
interstage  bleed  and  variable  inlet  guide  vanes.   An  after- 
burner section  was  present  but  was  not  used  during  the  test 
procedure.   The  engine  bleed  system  was  modified  to  make  the 
compressor  operate  with  less  than  normal  stall  margin.   The 
following  excerpt  is  quoted  from  reference  5: 

The  test  procedure  was  to  establish  an  engine 
operating  condition  in  the  vicinity  of  the  com- 
pressor stall  point.   The  FM  tape  recorder  was 
started  and  a  steady-state  data  scan  was  made. 
Then,  while  the  FM  tape  recorder  was  still  run- 
ning, the  compressor  pressure  ratio  was  increased 
slightly  and  the  throttle  adjusted  to  maintain  a 
constant  rotor  speed.   If  stall  did  not  occur, 
this  procedure  was  repeated  until  the  stall  event 
was  recorded  on  the  FM  tape.   In  some  cases  the 
stall  occurred  after  the  engine  had  stabilized 
at  an  equilibrium  operating  condition.   These 
were  called  'drift'  stalls. 
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All  eleven  of  the  stalls  examined  in  this  thesis  were  of 
the  "drift"  type.   NASA  Lewis  considered  only  seven  of  the 
eleven  cases  in  their  analysis.   Cases  6,  8,  9,  and  10  were 
excluded  on  the  grounds  either  that  their  causes  were  attri- 
butable to  steady-state  distortion  levels  or  that  they  showed 
very  little  loss  in  stall  margin.   Table  I  shows  steady-state 
operating  conditions  for  all  eleven  cases  presented  in  this 
thesis.   The  NASA  Lewis  analysis  identifies  each  case  number 
as  a  steady-state  reading  number.   These  NASA  identifying 
numbers  are  included  in  Table  I  for  comparison  [ref  5] . 

Steady-state  and  dynamic  pressure  instrumentation  at  the 
compressor  face  is  shown  in  figure  7.   Only  the  fluctuating 
component  of  pressure  was  recorded  from  each  of  the  dynamic 
probes.   Absolute  pressures  were  later  obtained  by  adding 
the  fluctuating  component  to  the  steady-state  value. 

Frequency  response  of  the  dynamic  probes  was  flat  to 
about  2000  hertz.   The  fluctuating  signal  from  each  pressure 
transducer  was  recorded  on  FM  multiplexed  tape  at  152.4  cen- 
timeters per  second.   The  data  in  a  225  millisecond  period 
just  prior  to  stall  were  digitized  at  an  effective  rate  of 
8000  points  per  second  per  channel.   These  data  when  processed 
by  the  computer  program  listed  in  the  back  of  this  thesis 
provided  the  distribution  of  total  pressure  across  the  com- 
pressor face  with  a  time  spacing  of  0.125  milliseconds  between 
distributions.   The  distributions  of  total  pressure  were  fur- 
ther processed  to  provide  distributions  of  both  circumferen- 
tial and  radial  vorticity  as  a  function  of  time.   In  order  to 
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TYPE    OF  PRESSURE 
o        STATIC 
o        TOTRL 

OPEN    SYMBOL      IMPLIES 
STEADY  STRTE 
SOLID   SYfABOL    IMPLIES 
DYNAMIC 


Z\Z.S 


Compressor   Pace 
FIGURE  7 

* 

reduce  the  amount  of  calculating  involved,  only  one  out  of 
five  readings  was  actually  processed  so  that  there  was  five 
times  0.125  or  0.625  milliseconds  of  time  between  readings 
actually  processed. 

The  determination  of  stall  time  was  made  by  NASA  Lewis. 
The  signals  from  the  compressor  static  pressure  transducers 
were  put  on  strip  charts  and  examined  by  NASA.   For  all 
cases,  the  row  of  instrumentation  at  the  290  degree  location 
showed  the  earliest  indication  of  stall.   Since  the  origin 
of  stall  in  a  compressor  stage  causes  an  abrupt  loss  in  the 
airflow  pumping  of  that  stage,  a  compression  wave  is 
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1   \ 


propagated  upstream.   This  compression  wave  is  referred  to 
as  a  hammershock.   The  following  excerpt  is  quoted  from  ref- 
erence 5: 

Because  there  were  only  two  axial  rows  of  in- 
strumentation inside  the  compressor,  it  was  not 
possible  to  determine  exactly  when  stall  originated 
inside  the  compressor.   Although  the  first  indica- 
tion was  seen  at  the  290  degree  location,  a  stall 
zone  could  originate  just  past  the  45  degree  loca- 
tion but  not  be  detected  until  it  rotated  past  the 
instrumentation  at  290  degrees.   Since  the  stall 
zone  rotated  at  about  half  rotor  speed  (8  msec/rev) 
there  could  be  as  much  as  a  6-millisecond  delay  be- 
fore it  was  detected  at  the  290  degree  location. 

A  one-millisecond  period  for  transport  time  between  sensors 

and  stall  stage,  plus  time  for  blade  response,  was  predicted 

by  NASA.   This  produces  a  six-plus-one  or  seven-millisecond 

possible  time  lag  between  actual  initialization  of  stall  and 

the  stall  time  listed  in  this  thesis. 
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V.   THREE  INDICES  APPLIED  TO  DATA 

A.   EXAMINATION  PROCEDURE 

The  computer  program  first  developed  by  Shoemaker  [ref 
3]  was  adapted  to  compute  each  index  as  a  function  of  time 
for  the  entire  time  interval  provided  by  NASA  Lewis.   This 
time  interval  varied  slightly  but  was  generally  about  215 
milliseconds  of  data  prior  to  first  indication  of  stall. 
The  DRAW  subroutine  of  the  Naval  Postgraduate  School  was  used 
to  plot  each  index.   These  plots  are  included  in  Appendix  A 
and  are  plots  of  the  instantaneous  version  of  each  index.   A 
numerical  printout  was  also  provided  by  the  program.   A  sam- 
ple numerical  printout  for  the  Shoemaker  index  case  one  is 
included  at  the  end  of  this  thesis.   Table  II  provides  an 
explanation  of  the  terms  for  the  numerical  printout. 

In  addition  to  the  computer  plots,  which  may  be  used  to 
determine  if  the  given  index  is  an  accurate  predictor  of 
stall,  a  simple  statistical  approach  was  used.   After  the 
first  fifty  values  of  the  index  were  computed,  the  mean  and 
standard  deviation  for  those  fifty  values  were  obtained. 
Then,  as  the  next  fifty  values  of  the  index  were  computed, 
the  following  test  was  made.   If  the  average  of  the  three 
previous  values  of  the  index  at  any  time  exceeded  a  certain 
criterion  value,  the  computer  program  printed  out  a  stall 
warning  with  the  time  in  milliseconds  of  occurrence.   The 
criterion  level  was  set  in  the  following  manner: 
Test  =  mean  +  2.5  *  standard-deviation. 
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N 

FACT(N)  = 
DFACT(N)  = 
TIME 

PTAVE 

MEAN 

STD-DEV   = 


TABLE  II 
EXPLANATION  OF  PROGRAM  TERMS 
index  number  (a  function  of  time), 
value  of  instantaneous  index  at  time  "N". 
average  of  three  previous  FACT(N)  values. 
"N"  time  in  milliseconds.   Starts  at  zero  for 
each  new  case. 

average  total  pressure  over  compressor  face  at 
time  "N". 

mean  value  of  FACT(N)  for  50  previous  values  of 
N,  when  not  zero. 

standard  deviation  of  FACT(N)  for  50  previous 
values  of  N,  when  not  zero. 


At  the  end  of  each  50  values  of  the  index,  a  new  mean  and 

standard-deviation  were  computed  and  the  test  criterion  was 

established  for  the  next  50  values  of  the  index.   This  simple 

test  was  used  on  all  thirty-three  runs  made,  and  the  results 

are  summarized  in  the  next  section.   Obviously,  the  test 

criterion  could  be  lowered,  with  resulting  greater  occurrence 

of  stall  prediction,  but  this  would  also  increase  the  number 

of  false  warnings.   The  value  of  test  was  initially  set  to 

predict  stall  for  the  K^,,  index  Case  three  without  any  false 
^  DM 

warnings.   The  test  criterion  could  have  been  set  to  an  opti- 
mum value  for  each  particular  index  with  possible  better  re- 
sults but  this  method  was  not  used.   The  same  test  was  applied 
to  all  cases.   The  statistical  method  lias  the  real  advantage 
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of  not  being  subject  to  interpretation.   All  thirty-three 
computer  plots  are  included  in  this  thesis,  since  they  cer- 
tainly are  subject  to  personal  interpretation  as  to  the  pre- 
diction of  compressor  stall. 

B.   COMPARISON  OF  RESULTS 

There  are  two  indicators  of  the  results  of  the  data  re- 
duction.  First  are  the  computer  plots,  which  provide  an 
eyeball  examination  of  results.   A  matrix  summation  of  the 
computer  plot  results  is  provided  by  Table  III.   The  second 
method  was  the  use  of  statistical  mean  and  standard-deviation 
discussed  in  the  preceding  section.   A  matrix  summation  of 
the  statistical  method  is  provided  by  Table  IV.   In  both 
tables,  where  stall  was  predicted,  the  time  to  stall  is  given 
in  milliseconds. 

There  are  several  points  of  interest.  As  had  been  pre- 
viously expected,  there  is  some  correlation  between  the  K  ., 
index  and  the  Downs  index,  since  both  rely  on  circumferential 
distortion.  Using  the  statistical  method,  all  of  the  cases 
predicted  by  the  K  „  index  were  also  predicted  by  the  Downs 
index.  The  Downs  index  predicted  two  additional  cases  uti- 
lizing this  method.   Using  the  computer  plot  interpretation, 

the  IC~W  method  predicts  six  cases  and  the  Downs  index  pre- 
DM 

diets  five,  with  four  predictions  common  to  both  indices. 

If  one  were  to  use  a  combination  of  both  vorticity  based 
indices,  the  statistical  method  would  have  predicted  eight 
out  of  the  eleven  cases  presented.   Using  this  method  there 
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TABLE  III 
COMPUTER  PLOT  SUMMARY 

(based  on  visual  interpretation) 

RUN     PRATT  &  WHITNEY   SHOEMAKER        DOWNS 
Koet   INDEX          INDEX            INDEX 

1 

No  Prediction 

22.5  millisecs 

23.0  millisecs 

2 

3.0  millisecs 
3  same  magnitude 

25.0  millisecs 

No  Prediction 
one  false  peak 

3 

3.0  millisecs 

21 .0  millisecs 
4  false  peaks 

6.0  millisecs 
2  false  peaks 

4 

4.0  millisecs 
2  larger  peaks 

No  Prediction 

No  Prediction 
one  false  peak 

5 

tit-   t>~,~  ;j  .s  „  x  ■<  ^v, 
iiv  x  j.  u u. x u  o luu 

^  «  U  ullxxlo&Co 

2  smaller  peaks 

No  Prediction 

6 

No  Prediction 

No  Prediction 

No  Prediction 

7 

No  Prediction 

No  Prediction 

No  Prediction 

8 

4.5  millisecs 

No  Prediction 

2.0  millisecs 

9 

4.0  millisecs 
one  false  peak 

16.0  millisecs 
5  same  magnitude 

2.0  millisecs 

10 

4.0  millisecs 

No  Prediction 

2.0  millisecs 

1 1 

No  Prediction 

6.0  millisecs 
2  larger  peaks 

No  Prediction 

totals 

6  Predictions 

6  Predictions 

__ _ 

5  Predictions 

32 


TABLE  IV 

STATISTICAL  SUMMARY 

RUN 

PRATT  &  WHITNEY 
K   INDEX 

SHOEMAKER 
INDEX 

DOWNS 
INDEX 

1 

No  Prediction 

21 .8  millisecs 

22.28  millisecs 

1  false  warning 

2  false  warnings 

2 

No  Prediction 

23.9  millisecs 

No  Prediction 
2  false  warnings 

3 

0.9  millisecs 

20.2  millisecs 

4.38  millisecs 

1  false  warning 

2  false  warnings 

4 

No  Prediction 

No  Prediction 

No  Prediction 

1  false  warning 

1  false  warning 

5 

No  Prediction 

No  Prediction 

No  Prediction 

6 

No  Prediction 

No  Prediction 

No  Prediction 

7 

No  Prediction 

No  Prediction 

19.33  millisecs 
1  false  warning 

8 

2.88  millisecs 

No  Prediction 
2  false  warnings 

2.33  millisecs 

9 

0.4  millisecs 

75.8  millisecs 
2  false  warnings 

0.0  millisecs 

10 

1 .62  millisecs 

No  Prediction 

2.33  millisecs 

1  false  warning 

1  false  warning 

11 

No  Prediction 

11.8  millisecs 

No  Prediction 

totals 

4  Predictions 

5  Predictions 

6  Predictions 
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would  have  been  seventeen  false  warnings.   The  number  of 
false  warnings  could  be  reduced  by  a  reduction  in  the  test 
criterion  to  a  point  where  no  loss  in  stall  prediction  would 
occur. 

There  could  be  as  much  as  seven  milliseconds  of  time  lag 
between  actual  stall  time  and  predicted  stall  time.   This 
presents  a  problem  since  the  data  could  include  the  effects 
of  the  hammershock  which  occurs  after  stall,  if  the  time  lag 
is  too  great.   It  is  certainly  no  problem  to  produce  an  index 
which  will  sense  the  hammershock  and  the  question  arises,  did 
the  index  actually  predict  stall  or  were  the  effects  of  the 
hammershock  being  felt?   By  examining  the  numerical  printout 
of  average  total  pressure  as  a  function  of  time  and  taking 
into  account  the  possible  seven-millisecond  lag  time,  it  was 
apparent  that  this  could  have  been  a  problem  in  only  four 
cases.   These  four  cases  should  be  held  suspect  and  are  cases 
3,  8,  9,  and  10.   It  is  noteworthy  that  the  Shoemaker  index 
predicted  stall  in  two  of  these  four  cases  sufficiently  far 
in  advance  to  preclude  the  problem. 
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VI.   CONCLUSIONS 

An  investigation  into  the  problem  of  predicting  compres- 
sor stall  using  instantaneous  distortion  parameters  was  made. 
Some  success  was  evidenced  by  both  the  statistical  method 
and  computer  plot  interpretation,  but  it  was  not  possible  to 
predict  stall  in  every  case  using  one  or  any  combination  of 
the  indices  explored.   This  seems  to  indicate  that  an  axial 
compressor  is  subject  to  many  instantaneous  distortion  pat- 
terns as  well  as  to  steady-state  distortion,  whereas  one  in- 
dex may  only  search  for  one  category  of  distortion.   If  the 
cause  of  compressor  stall  was  in  every  case  instantaneous 
radial  distortion  over  the  outer  portion  of  the  compressor 
face,  then  in  every  case  the  Shoemaker  index  should  predict 
stall.   If  the  cause  was  in  every  case  circumferential  dis- 
trotion  over  some  arbitrary  sixty  degree  sector  then  the 
Downs  index  would  be  successful.   If  the  distortion  conforms 
to  the  ring  pattern  explored  by  the  Pratt  and  Whitney  index, 
then  that  method  would  succeed.   In  fact,  it  seems  naive  to 
believe  that  compressor  stall  could  be  predicted  in  every 
case  using  any  one  index.   The  indices  explored  in  this  the- 
sis and  the  indices  used  by  industry  today  may  certainly  be 
improved  upon,  but  it  is  doubtful  that  any  one  index  will 
prove  to  be  a  panacea. 

In  spite  of  the  compound  nature  of  compressor  stall,  in- 
stantaneous indices  are  useful.  The  Pratt  and  Whitney  index 
used  in  this  thesis  was  developed  for  the  F-111A  with  a 
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specific  inlet-engine  combination.   It  is  in  that  context 
that  distortion  parameters  would  appear  most  useful  -  that 
is,  to  examine  an  inlet-engine  system  for  compatibility. 
The  performance  of  an  inlet-engine  system  which  too  often 
experiences  one  type  of  distortion  pattern  may  be  improved 
upon;  distortion  indices  would  be  very  useful  in  isolating  a 
predominant  pattern. 

The  use  of  vorticity  as  an  index  parameter  has  some  real 
advantages.   Both  temperature  and  pressure  gradients  may  be 
integrated  into  a  single  parameter.   Vorticity  may  be  visual- 
ized more  simply  and  is  a  natural  characteristic  of  the  fluid 
The  Shoemaker  index,  when  successful,  harbingers  stall  suffi- 
ciently far  in  advance  possibly  to  enable  correction  of  the 
deficiency  using  engine  controls.   Vorticity  appears  to  be  a 
promising  dimensionless  parameter  and  deserves  further  inves- 
tigation . 
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APPENDIX  A 
COMPUTER  PLOTS 

The  following  thirty-three  computer  plots  provide  the 
instantaneous  version  of  each  of  the  three  indices  considered 
In  all  cases,  the  horizontal  scale  is  read  in  milliseconds. 
The  vertic.nl  scale  gives  the  index  value  and  is  a  non-dimen- 
sional number.   No  consideration  was  given  to  index  magnitude 
in  predicting  stall.   For  any  given  case,  a  strong  peak  just 
before  stall  in  relation  to  that  particular  case  was  consid- 
ered to  be  sufficient  for  stall  prediction. 
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